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OBJECTIVE — Two meta-analyses of randomized trials of statins found increased risk of type 2 
diabetes. One possible explanation is bias due to differential survival when patients who are at 
higher risk of diabetes survive longer under statin treatment. 

RESEARCH DESIGN AND METHODS— We used electronic medical records from 500 
general practices in the U.K. and included data from 285,864 men and women aged 50-84 years 
from January 2000 to December 2010. We emulated the design and analysis of a hypothetical 
randomized trial of statins, estimated the observational analog of the intention-to-treat effect, and 
adjusted for differential survival bias using inverse-probability weighting. 

RESULTS — During 1.2 million person-years of follow-up, there were 13,455 cases of type 2 
diabetes and 8,932 deaths. Statin initiation was associated with increased risk of type 2 diabetes. 
The hazard ratio (95% CI) of diabetes was 1.45 (1.39-1.50) before adjusting for potential con- 
founders and 1.14 (1.10-1.19) after adjustment. Adjusting for differential survival did not 
change the estimates. Initiating atorvastatin and simvastatin was associated with increased risk 
of type 2 diabetes. 

CONCLUSIONS — In this sample of the general population, statin therapy was associated 
with 14% increased risk of type 2 diabetes. Differential survival did not explain this increased 
risk. 
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Recent analyses of randomized clin- 
ical trials of statins have shown a 
higher incidence of diabetes in the 
treatment arm (1,2). A meta-analysis of 13 
placebo-controlled randomized trials 
showed 9% (95% CI 2-17%) higher odds 
of diabetes in the statin arm. Another meta- 
analysis of five secondary prevention trials 
that compared more- versus less-intensive 
statin treatment reported 12% (4-22%) 
higher odds of diabetes in the more-intensive 
treatment arms (2). One observational 
study also reported a positive association 
between statin use and the subsequent 
risk of diabetes (3), although another study 
reported different associations (null, posi- 
tive, or negative) depending on the par- 
ticular statins used (4). The collective 



evidence led the Food and Drug Adminis- 
tration to require that the labels of all statin 
drugs be changed to indicate the potential 
increased risk of diabetes as a side effect. 

Part of the observed association be- 
tween statin use and risk of diabetes in both 
randomized trials and observational stud- 
ies might be attributable to bias from 
differential survival. Because statins reduce 
the risk of death from cardiovascular dis- 
eases (CVD) (5), survivors who were treat- 
ed with statins could have a worse risk 
profile for CVD than untreated survivors. 
Therefore, the treated survivors might be at 
higher risk of diabetes simply because CVD 
risk factors such as overweight/obesity, 
smoking, and poor diet are also risk fac- 
tors for type 2 diabetes (Fig. 1). Indeed, a 
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meta-regression analysis of randomized tri- 
als found a stronger harmful effect in trials 
that enrolled older participants or those 
with a higher baseline mean BMI, which 
suggests potential bias due to differential 
survival (1). Neither the randomized trials 
nor the observational studies addressed this 
possibility, perhaps because the absolute 
difference in survival in the randomized tri- 
als was only 1.4 percentage points (1). 

We set out to examine the role of bias 
due to differential survival by applying 
inverse probability weights to data from a 
database of electronic medical records 
from the U.K. 

RESEARCH DESIGN AND 
METHODS— We emulated the design 
and analysis of a randomized trial of statins 
using observational data as described in 
detail elsewhere (6). We used data from the 
Health Improvement Network (THIN), 
which is a large computerized database of 
anonymized longitudinal medical records 
from >500 primary care practices in the 
U.K. The database includes demographic 
information as well as symptoms and diag- 
noses, summaries of referrals and hospital- 
izations, outpatient prescriptions, test 
results, and lifestyle information for every 
registered patient. 

We started the follow-up on 1 January 
2000 and followed patients through 31 
December 2010. We included men and 
women 50-84 years of age who did not 
have diabetes (type 1, type 2, or gestational 
diabetes), had not been prescribed any sta- 
tins in the past 2 years, had at least 2 years 
of continuous recording in the database 
(operationalized as at least 2 years since 
the first recorded prescription for any 
drug), and at least one health contact 
within the past 2 years. We excluded pa- 
tients who had cancer, chronic liver, or kid- 
ney disease, schizophrenia, or used 
antipsychotic drugs at baseline. 

The intervention of interest was 
being prescribed any statin (atorvastatin, 
fluvastatin, pravastatin, rosuvastatin, and 
simvastatin). As several previous studies 
had indicated that effects may depend on 
the type of statins (7), we also compared 
the effects of different statins versus no 
treatment. 
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Figure 1 — Simplified directed acyclic graph 
presenting the potential for differential sur- 
vival Di, death during follow-up (e.g., due to 
cardiovascular disease); DM 2 , diabetes melli- 
tus by the end of follow-up; L 1} a vector of 
common risk factors of diabetes and mortality 
such as age, obesity, and smoking; S 0 , statin 
therapy at baseline. 

The outcome of interest was incident 
type 2 diabetes. The diagnosis was based on 
the Read classification codes assigned by 
the general practitioners or use of hypogly- 
cemic drugs or insulin; the date of incident 
diabetes was the earliest of diabetes codes 
or diabetes medication. For the majority of 
cases, the type of diabetes was specifically 
reported by the physician. If the physician 
used an unspecific diagnostic code (e.g., 
diabetes mellitus), we reviewed the pa- 
tient's medical record back to 1 year before 
the diagnosis including any referral letters 
and physicians' free-text comments to as- 
sign the type of diabetes. If the age of onset 
was <35 years and the patient had 1 or 
more prescriptions for insulin and <1 
year of oral hypoglycemic treatment, the 
case was classified as type 1 diabetes. Con- 
versely, if the age of onset was >50 years 
and the patient used oral hypoglycemic 
treatment for at least 1 year, the case 
was classified as type 2 diabetes (8-10). 
A previous THIN study with a similar di- 
abetes ascertainment algorithm estimated a 
diabetes prevalence that closely matched 
the prevalence in the Health Survey of 
England, which is a national population 
survey (9). 

We also examined the effect of statins 
on peptic ulcer as a negative outcome 
control (11) because there is no plausible 
mechanism that statins would affect risk 
of peptic ulcer, and all of the observed 
association between statin use and inci- 
dent peptic ulcer should be due to bias. 
A validation study of uncomplicated pep- 
tic ulcer cases in THIN found a positive 
predictive value of 94% (12). 

Potential confounders included age, 
sex, Townsend deprivation index (a mea- 
sure of neighborhood deprivation), tobacco 
smoking (never, former, or current), alcohol 
use, test results (LDL and HDL cholesterol, 
serum triglycerides), physical examination 
results [BMI, systolic blood pressure (SBP)], 
history of diabetes in first-degree relatives, 



concurrent chronic diseases, use of other 
drugs, and number of visits, referrals, and 
hospitalizations in the past 3 months. See 
Table 1 for the complete list of variables. We 
truncated the distribution of risk factors 
at a threshold based on expert opinion 
to prevent implausible values (0.2%) from 
affecting the results. We carried the last ob- 
servation forward for a maximum of 12 
months for LDL, HDL, and SBP and of 24 
months for alcohol use, smoking, and BMI. 

Statistical analyses 

The details of the methods are described 
elsewhere (6). Briefly, we generated a set 
of emulated trials in which the first trial 



recruited participants for the duration of 
the first month of the study (in this study, 
January 2000). Eligibility criteria were ap- 
plied, and patients who were eligible were 
categorized as either initiators (i.e., those 
who started using statins in January 2000) 
or noninitiators (i.e., those who did not 
start statin therapy in January 2000). 
Noninitiators were eligible for the emu- 
lated trial recruiting in February 2000, 
whereas initiators should wait for at least 
24 months after they discontinued statin 
therapy to become eligible again. We con- 
structed 131 emulated trials using THIN, 
pooled the data from all trials, and 
followed patients until occurrence of 



Table 1 — Characteristics of initiators and noninitiators of statin therapy at the start of the 
triaVs follow-up: THIN trials 2000-2010 





Noninitiators 






(3,765,906 


Initiators (48,013 


Characteristic 


person-trials) 


person-trials) 


Male (%) 


44.2 


50.2 


Townsend deprivation score of >3 (%) 


44.2 


48.0 


Age (years) 


61.3(9.3) 


64.1 (8.4) 


LDL cholesterol (mmol/L) 


3.48 (0.86) 


4.18(0.95) 


HDL cholesterol (mmol/L) 


1.51 (0.43) 


1.44 (0.41) 


BMI (kg/m 2 ) 


27.7(5.1) 


28.1 (4.9) 


SBP (mmHg) 


137(16) 


143 (18) 


Alcohol use (units/day) 


0.06 (0.24) 


0.06 (0.24) 


Doctor visits in the past 3 months 


3.8 (7.3) 


5.6 (5.2) 


Referrals in the past 3 months 


(J./ u./j 


l.U u./j 


Hospitalizations in the past 3 months 


0.04 (0.3) 


0.06 (0.2) 


Alcoholism (%) 


1.9 


2.1 


Smoking prevalence (current or former) (%) 


41.7 


49.9 


Hypertension (%) 


50.7 


62.1 


Antihypertensive use (%) 


47.9 


66.8 


NSAIDs use (%) 


9.8 


12.7 


Aspirin use (%) 


7.8 


26.7 


Other lipid-lowering drugs use (%) 


0.7 


0.8 


(3-Blockers use (%) 


14.8 


22.8 


Hormone replacement therapy (% of women) 


9.8 


5.8 


Chronic obstructive pulmonary disease (%) 


3.4 


4.7 


Oral steroids use (%) 


1.8 


2.0 


Inhaled steroids use (%) 


5.7 


6.5 


Atrial fibrillation (%) 


1.6 


2.2 


Depression (%) 


21.0 


21.2 


Antidepressant use (%) 


8.4 


9.6 


Hypothyroidism (%) 


7.4 


7.6 


Osteoporosis (%) 


3.8 


3.9 


History of transplant (%) 


0.03 


0.03 


Immunosuppression therapy (%) 


0.4 


0.4 


Psoriasis (%) 


4.0 


4.2 


Rheumatoid arthritis (%) 


2.0 


2.2 


Chemotherapy (%) 


0.8 


0.8 


Radiotherapy (%) 


0.04 


0.06 


Chronic pancreatitis (%) 


0.05 


0.07 



Numbers are mean (SD), unless otherwise specified. NSAID, nonsteroidal anti-inflammatory drug. 



care, diabetesj ournals . org 



Diabetes Care, volume 36, May 2013 



1237 



Statins and diabetes, evaluating survival bias 



diabetes, death, or administrative end of 
follow-up, whichever occurred earlier. 

We estimated the observational analog 
of the intention-to-treat hazard ratio (HR) 
of diabetes for initiators versus noninitia- 
tors by fitting a Cox model that included 
the baseline confounders for each trial. We 
used a robust variance estimator because 
each individual may participate in multiple 
trials. 

To evaluate the sensitivity of our esti- 
mates to differential survival, we used in- 
verse probability weighting to account for 
censoring by death. In each trial, we esti- 
mated the monthly probability of death 
from any cause using a logistic regression 
conditional on common risk factors of 
death and diabetes. We carried forward 
the last observed value of the selected 
covariates (listed in Table 1) indefinitely 
and added a set of additional variables 
that indicated the length of time passed 
since last measurement of SBP and LDL 
and HDL cholesterol levels. The denomi- 
nator of the stabilized weight for month m 
was calculated as the estimated probability 
of surviving from baseline until month m 
(6). We truncated the distribution of 
weights at 50 (which is larger than the 
99th percentile) to prevent the influence 
of observations with very large weights (6). 

Informally, the weights create a pseu- 
dopopulation in which censoring by death 
is independent of the measured risk factors. 
Although such pseudopopulation does not 
exist in practice, this approach assesses the 
sensitivity of the effect estimates to differ- 
ences in risk factor profiles between the 
surviving treated and untreated. If the 
weighted analysis yields similar estimates 
for the effect of statins on diabetes risk as 
the unweighted analysis, strong bias due to 
differential survival is unlikely. We hypoth- 
esized that any bias due to differential 
survival would be larger in the older sub- 
groups and conducted an analysis re- 
stricted to patients >62 years of age (the 
median age at baseline). We also examined 
the effect of each statin drug compared with 
no statins and the effect of hydrophilic sta- 
tins (pravastatin, rosuvastatin, or fluvasta- 
tin) and lipophilic statins (simvastatin or 
atorvastatin) compared with no statins sep- 
arately. 

RESULTS — Figure 2 shows a flow chart 
of the process for selecting eligible study 
participants, and Table 1 presents their 
characteristics by statin initiation status. 
We included data from 285,864 patients 
who did not have diabetes at baseline. On 
average, the 48,013 statin initiators were 



2.8 years older, had 0.4 kg/m 2 higher BMI, 
and 0.7 mmol/L higher LDL than the 
3,765,906 noninitiators. A larger propor- 
tion of statin initiators used antihyperten- 
sive medications and aspirin, and they had 
on average two more doctor visits in the past 
3 months compared with noninitiators. 

During 1.17 million person-years of 
follow-up, there were 13,455 new cases of 
diabetes and 8,932 deaths. Median time 
of follow-up was 28 months among initiators 
and 29 months among noninitiators. The 
crude incidence rate of type 2 diabetes was 
15.9 per thousand person-years among 
initiators and 11.3 among noninitiators. 
The results of the intention-to-treat analysis 
are presented in Table 2. The unadjusted HR 
(95% CI) was 1.45 (1.39-1.50). Adjustment 
for confounding by indication reduced the 
HR to 1.14 (1.10-1.19). The adjusted HR 
was 1.24 (1.14-1.34) in the first year of 
follow-up and 1.12(1.06-1.17) after exclud- 
ing the first year. 

By the end of follow-up, there were 78 
unique deaths among 14,019 statin initia- 
tors and 8,854 among 271,845 noninitia- 
tors. Initiating statin therapy reduced the 
risk of death significantly. The adjusted 
intention-to-treat HR for mortality was 
0.88 (0.84-0.93). The estimated inverse- 
probability weights had a mean (SD) of 
1.00 (0.23). The HRs of diabetes in the 
weighted analyses were almost identical 
to the unweighted ones both in the main 



analysis and in the analysis restricted to 
older patients. 

Of all statin initiators, 79.0% used 
simvastatin and 15.5% atorvastatin. Table 
3 presents the adjusted and weighted HR 
of diabetes for initiation of each statin or a 
group of statins versus no statin therapy. 
Lipophilic statins (i.e., simvastatin and 
atorvastatin) were associated with an in- 
creased risk of diabetes. The HR was 1.14 
(1.09-1.20) for simvastatin and 1.22 
(1.12-1.32) for atorvastatin. There were 
not enough events to conduct a meaning- 
ful comparison of initiation of other sta- 
tins. The intention-to-treat analysis 
among those >62 years of age gave an 
HR of 1.16 (1.10-1.23). Finally, the HR 
of peptic ulcer comparing statin initiation 
versus no initiation was 0.98 (0.91-1.06). 

CONCLUSIONS— We estimated that 
statin initiators had a 14% increased in- 
cidence of diabetes compared with non- 
initiators in this general population in 
the U.K. The increased risk was observed 
with simvastatin and atorvastatin, the 
most commonly used statins in this pop- 
ulation. We found no evidence that the 
increased risk may be explained by differ- 
ential survival of people treated with 
statins. Our estimates did not change after 
inverse-probability weighting for death or 
after restricting the analysis to older indi- 
viduals. 



8.9 million potential 
person-trials 
(1,192,384) 



3,813,919 eligible 
person-trials 
(285,864) 



-16.2 million major chronic diseases 
- 66. 1 million no recent data on 
confounders 

-2.8 million prior statin use or <2 
years on database 



48,013 initiators 
(14,019) 



- 1658 died (78) 

- 4632 were lost to follow-up 



3,765,906 non-initiators 
(271,845) 



- 2890 diabetes cases (902) 

- 38,833 alive and event-free at 
the end of follow-up 



- 100,739 died (8854) 

- 334,709 were lost to follow-up 



- 149,746 diabetes cases (12,543) 

- 3,180,712 alive and event-free at 
the end of follow-up 



Figure 2 — Flow chart of person-trials in the analysis: THIN database 2000-2010. Numbers in 
parentheses indicate unique individuals. 
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Table 2— HR (95% CI) for the analog of the 
intention-to-treat effect of statins on type 2 
diabetes: THIN database 2000-2010 





ITT, HR 


Analysis 


(95% CI) 


Noninitiators 


Reference 




Unadjusted 


1.45 (1.39-1.50) 


Unweighted, adjusted 




for baseline variables 3 


1.14(1.10-1.19) 


Weighted, 13 adjusted 




for baseline variables 3 


1.14(1.09-1.19) 



There were 13,455 unique cases of diabetes after 
1.17 million person-years of follow-up. ITT, in- 
tention-to-treat. a Adjusted for baseline covariates in 
Table 1. b The models used to estimate the denom- 
inator of the weights included the same variables in 
Table 1 during the follow-up as well as time since 
measurement of blood pressure and LDL and HDL 
cholesterol. 

Our study is observational, and there- 
fore, our estimates may be affected by 
unmeasured confounding. However, the 
magnitude of our estimate of the effect of 
statins on diabetes and total mortality are 
similar to those reported in meta-analyses of 
randomized trials (1,13). Also, our results 
cannot be fully explained by confounding 
by undiagnosed diabetes as the increased 
risk was detected even after the first year 
of follow-up. Our estimates might also be 
biased if the association between statins and 
diabetes were explained by better monitor- 
ing of statin users. However, we found the 
expected null association between statin ini- 
tiation and peptic ulcer, another condition 
that would be more likely to be diagnosed 
among better-monitored patients. 



Imperfect adherence to statin therapy 
could have attenuated our intention-to- 
treat effect estimate. However, the HR of 
diabetes in an as-treated analysis (6) for 26 
months of statin use, which was the average 
duration of treatment among initiators, ver- 
sus no statin use was 1.15 (results not 
shown). 

Our inverse-probability weighted esti- 
mates suggest that the differential survival 
of statin users does not explain their in- 
creased risk of diabetes. However, we may 
not have been able to adjust for all common 
risk factors of mortality and diabetes. For 
example, the THIN database does not in- 
clude any information on diet or physical 
activity. Furthermore, the data on smok- 
ing and alcohol use are not frequently 
updated which will result in measurement 
error and residual confounding by these 
factors. 

There are several proposed biological 
mechanisms for the harmful effects of statins 
on diabetes. Simvastatin and atorvastatin 
may increase insulin resistance (14-16), 
although fluvastatin may decrease it (17). 
Other potential pathways are through re- 
duction in synthesis of various metabolites 
such as isoprenoid, which leads to lower 
glucose uptake in adipocytes or reduced 
insulin secretion by blocking the calcium 
channels in the (3-cells (18). Further re- 
search is warranted to identify the role of 
various molecular pathways. 

It is worth noting that even with the 
slight increase in the risk of diabetes 
among statin initiators, the cardiovascular 
and mortality risk reduction outweighs 
this additional risk in patients with mod- 
erate or high CVD risk and even those 
who are at higher risk of diabetes (19). 



Table 3 — The number (proportion) of initiators for each statin or group of statins and 
adjusted and weighted intention-to-treat HRs (95% CI) for type 2 diabetes: THIN database 
2000-2010 



Statin used a 


Initiators 13 (%) 


Unique cases 
among initiators 


ITT, HR 
(95% CI) 


Simvastatin 


37,915 (79) 


584 


1.14(1.09-1.20) 


Atorvastatin 


7,438 (15.5) 


232 


1.22 (1.12-1.32) 


Pravastatin 


1,442 (3.0) 


45 


1.01 (0.84-1.21) 


Rosuvastatin 


1,175 (2.5) 


28 


1.11 (0.89-1.38) 


Fluvastatin 


282 (0.6) 


15 


1.02 (0.69-1.50) 


Lipophilic statins (simvastatin 
or atorvastatin) 


45,218 (94) 


814 


1.16 (1.11-1.21) 


Hydrophilic statins (pravastatin, 
rosuvastatin, or fluvastatin) 


2,895 (6) 


88 


1.04 (0.91-1.19) 



ITT, intention-to-treat. a Lovastatin is not available in the U.K. b A total of 239 out of 48,013 initiators used 
more than one statin at initiation time and contributed to more than one of these analyses. 



In summary, our results support a 
slightly increased risk of diabetes due to 
statin therapy in this general population 
in the U.K. We did not find any evidence 
of bias due to differential survival. Similar 
analyses could be conducted on the ran- 
domized trials to examine the possibility 
of this type of bias. 
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